function [Vmax]=question_2(alpha,h,T_set_2,range)
%问题二求解
    s=0;i=1;N=1000;
    while(s<N)
        v=65+35*rand;

        v=v/60;
        x=0:0.1:range(end);
        q2_time=x/v;
        q2_time=q2_time';

        n=length(q2_time);

        %设置问题二的环境温度
        T0=zeros(size(q2_time));
        for i=1:n
            T0(i)=get_K(T_set_2,x(i));
        end

        Divide_line=zeros(6,1);%存储各个温区的横坐标位置
        for i=1:4
            m=find(T0==T_set_2(i));
            Divide_line(i+1)=m(end);
        end
        Divide_line(6)=n;

        u0=25*ones(3,1);
        mid_T=zeros(size(q2_time));
        for i=1:5
            [mid_T(Divide_line(i)+1:Divide_line(i+1)),u0,u(:,Divide_line(i)+1:Divide_line(i+1)),~]=PDE(alpha(i),h(i),u0,T0(Divide_line(i)+1:Divide_line(i+1)),q2_time(Divide_line(i)+1:Divide_line(i+1)));
        end


        [flag]= limit(q2_time,mid_T);

        
        if(flag==1)
            V(i)=v;
            Vmax(i)=max(V);
            if(i==1)
                v0=V(1);
            else
                if(Vmax(i)>Vmax(i-1))
                    v0=v;s=0;
                else
                    s=s+1;
                end
            end
            i=i+1;
        end
    end
    
    figure
    plot(1:length(Vmax),Vmax,'k');
    xlabel('迭代次数')
    ylabel('最大速度/cm/min')
    title(['v_{max}=',num2str(v0)])

    figure
    plot(t,T,'r')
    xlabel('时间/s')
    ylabel('温度/℃')
end

